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Abstract 



Aims. We studied the thermal properties and chemical composition of the X-ray emitting plasma of a sample of bright members of the Taurus 
Molecular Cloud to investigate possible differences among classical and weak-lined T Tauri stars and possible dependences of the abundances 
on the stellar activity level and/or on the presence of accretion/circumstellar material. 

Methods. We used medium-resolution X-ray spectra obtained with the sensitive EPIC/PN camera in order to analyse the possible sample. The 
PN spectra of 20 bright (Lx ~ 10^" - 10^^' erg s"') Taurus members, with at least ~ 4500 counts, were fitted using thermal models of optically 
thin plasma with two components and variable abundances of O, Ne, Mg, Si, S, Ar, Ca, and Fe. Extensive preliminary investigations were 
employed to study the performances of the PN detectors regarding abundance determinations, and finally to check the results of the fittings. 
Results. We found that the observed X-ray emission of the studied stars can be attributed to coronal plasma having similar thermal properties 
and chemical composition both in the classical and in the weak-lined T Tauri stars. The results of the fittings did not show evidence for 
correlations of the abundance patterns with activity or accretion/disk presence. The iron abundance of these active stars is significantly lower 
than (~ 0.2 of) the solar photospheric value. An indication of slightly different coronal properties in stars with difi^erent spectral type is found 
from this study. G-type and early K-type stars have, on average, slightly higher Fe abundances (Fe ~ 0.24 solar) with respect to stars with later 
spectral type (Fe ~ 0.15 solar), confirming previous findings from high-resolution X-ray spectroscopy; stars of the former group are also found 
to have, on average, hotter coronae. 
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1. Introduction 

Studies on the chemical composition of stellar atmospheres are 
important to address several issues in astrophysics, such as the 
evolution of stellar galactic populations, the enrichment of the 
interstellar medium, the influence of elemental abundances on 
star formation processes and on the structure of stellar interiors 
through opacity effects. In the framework of coronal physics, 
solar and stellar observations in the UV and X-ray bands dur- 
ing the last two decades have suggested abundance anomalies 
that, if understood, could provide powerful new diagnostics for 
the physics of stellar out er atmospheres. In particular, previous 
observations of the Sun ( Meveij|l985t iFeldmanll 1 992 ). showed 
different chemical abundances in corona with respect to the 
photosphere, with a pattern depending on the First Ionization 
Potential (FIP) of the various elements. The present view is 
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that elements with FIP < 10 eV (low-FIP elements Mg, Si, Ca, 
Ni, Fe) are overabundant in the solar corona with respect to the 
photosphere by a factor ~ 4 - 5; elements with FIP > lOeV 
{high-FIP elements C, S, O, Ne, Ar, with C and S lying close 
to the boundary of this definition) have coronal abundances 
similar to the corresponding photospheric values. This solar- 
type FIP effect is most clearly present in long-lived active re- 
gions, whi le it is much less marked (or even absen t) in coro- 
nal holes dWiding & FeldmanI IT992I: ISheeIevl ll996'). In order 
to compare solar an d spatially-integrated stellar observations. 
Laming et alj ( 119941) investigated full-disk solar EUV spectra, 
and detected the FIP effect also in this case. 

Analyses of EUV line emission spectra from the coro- 
nae of low- and intermediate-activity stars (o-Centauri, ^ 
Bootis, e Eridani, observed with the EUVE satellite) revealed 
solar-like abundance anomalies, suggestin g that the FIP ef- 
fect is not an isolated solar phenomenon (Laming et al. I ll996l: 
Drake et alJ[l997l:lLaming & Drakelll999h . New studies, based 
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on XMM-Newton and Chandra high-resolution X-ray spectra, 
have led to claims of different FlP-dependent abundance pat- 
terns in the coronae of other stars: in particular, observations 
of very active stars, both single stars and active binary sys- 
tems, suggested a dependence of coronal abundances on FIP 
which is reversed with respect to the Sun and the other low- 
activity stars mentioned above, i.e. a depletion of low-FIP el- 
ements with respect to high-FIP elements (inverse FIP effect; 
terinkman et al."200lVAudard et al.||2003i iRobrade & SchmittI 
[2005)). Telleschi et al. (2005) derived coronal abundances for 
6 solar analogs (i.e. G-type main sequence stars) with different 
ages and rotation periods, and hence activity levels: this work 
suggested a change from an inverse FIP effect to a solar-type 
FIP effect with increasing age, and hence with decreasing ac- 
tivity. 

On the contrary, other authors found evidence for nei- 
ther a FIP nor an inverse FIP effect in t he coronae of active 
late-type stars ("Sa nz-Forcada et aklEoOsl: lArgiroffi et~aDl2004 ; 



Scelsi et al. 2005), but instead a minimum in the abundance 
pattern occurring at FIP ~ 10 eV, with increasing abundances at 
both lower and higher FIP values. In any case, a systematic re- 
sult emerging from stellar abundance studies is that the coronae 
of high-activity stars appear significantly metal -poor relative to 
the solar photosphere, with typical ir on abundances -0.1 -0.3 
times the lAnders & Grevesse lll989h solar photospheric value. 

A major problem encountered in coronal abundance stud- 
ies is related to data analysis issues. The requirements for an 
accurate analysis are good quality X-ray spectra (i.e. high S /N 
ratios), a careful selection of emission lines (not affected by 
blends and independent from the plasma density), accurate line 
emissivity functions describing the physics of atomic transi- 
tions, and a precise knowledge of the instrumental response. 
All these difficulties introduce errors that are often not ac- 
counted for in the reported (just statistical) measurement un- 
certainties, which therefore are likely underestimated in many 
studies. 

While the coronal metal deficiency of active stars seems 
to be a rather solid result, the question whether coronal abun- 
dances are different from the photospheric values, with a pat- 
tern possibly dependent on the FIP, is still debated and de- 
serves further investigation, especially because photospheric 
abundances of active, rapid rotators are mostly unavailable. 
Better observational evidence is definitely required to test the- 
oretical models proposed to explain the physical mechanism(s) 
responsible for such differences (e.g. ISchwadron et al.lll999t 
[Laming 2004). 

In classical T Tauri stars (CTTSs) material from the cir- 
cumstellar disk accretes onto the stellar surface and shock 
waves at the base of the accretion funnel are expected to heat 
the gas to X-ray emitting temperatures of a few million de- 
grees. From recent analyses of XMM-Newton and Chandra 
spectra of pre-main sequence stars, interesting peculiarities in 
the abundances of some CTTSs were also noted: TW Hydrae, 
TWA 5, and BP Tauri show anomalous Ne/Fe abundance ratios 
(-5-10 solar. IStelzer & SchmittlEooi: lArgiroffi et"aD l2005: 



formed, but they might also be due to some process of chemical 
segregation during the evolution of circumstellar disks. In fact, 
most of the soft X-ray emission of TW Hydrae and part of the 
emission of BP Tauri have been attributed to accretion shocks 
(while the situation for TWA 5 is more uncertain). Iron and 
other metals in the disk could conde nse in grains , as predicted 
by so me model s of cloud chemistry (ISavage & Sembach 199S 
[Charnlevl[l997i) . If we suppose that the gas and dust phases in 
the circumstellar disk are separated for some reason, with the 
metals remaining at larger distance from the star, while neon 
and other noble elements in the gas phase accrete more effi- 
ciently onto the stellar surface, then the gas heated by the shock 
at X-ray emitting temperatures is expected to be metal-deficient 
and enriched in neon. The high Ne/O ratio of TW Hydrae may 
be another sign of such a segregation pro cess occurred in th e 
circumstellar disk of this star, as argued bv lDrake et al.l ( l2005h . 
Should this scenario be proven correct for the CTTSs and the 
gas-dust separation mechanism understood, the chemical com- 
position of the hot plasma could provide new insight into the 
study of circumstellar disks. The relevance for our understand- 
ing of planetary formation could also be significant. 



Most recently. lMaggio et al.l (l2007h studied a sample of 146 
X-ray bright sources associated with members of the young 
Orion Nebula Cluster (ONC), and found that their emission 
originates from a plasma with temperatures and elemental 
abundances very similar to those of active coronae in older 
stars. Moreover, the abundance distributions are compatible 
with a single pattern vs. FIP for all stars for all elements 
with the only possible exception of the Calcium abundance. 
Comparison with abundances in stellar photospheres and in the 
gaseous component of the nebula indicate a significant under- 
abundance at coronal level only for iron (by a factor 1.5-3), 
while there is agreement (and hence no evidence of a FIP- 
related effect) for all the other elements with available mea- 
surements. 

In another recent study based on nine high-resolutio n RGS 
spectra from the XEST survey, ' Telleschi et al. found 
similar abundances in accreting and non-accreting stars and in- 
dication of higher Fe abundances in the G-type stars of their 
sample compared to stars with later spectral type. They also 
found that CTTSs display a soft excess when compared to 
WTTSs or zero-age main-sequence stars, and that high electron 
densities and Ne/Fe abundance ratios are not characteristic of 
all accreting PMS stars. 

In this paper we analyze X-ray spectra taken with XMM- 
Newton to derive the thermal properties and chemical compo- 
sition of the X-ray emitting plasma of a sample of 20 young 
stellar objects (7 CTTSs and 13 WTTSs) belonging to the 
Taurus-Auriga star forming region and spanning an age range 
of ~ 1 - 10 Myr. This region has been recently the s ubject of the 



XMM -Newton Extended Survey of Taurus (XEST, iGiidel et al 



Robrade & Schmit t 2006), significantly higher than those typi- 
cally found for other stars (A(Ne)/A(Fe) ~ 3). Such peculiari- 
ties could be inherited from the original cloud where these stars 



[20071) . which allowed us to study the X-ray emission from more 
than 150 classical and weak-lined T Tauri stars in the Taurus 
Molecular Cloud (TMC). The present work focuses on some 
of the most active members of the TMC, with the aim to under- 
stand whether accreting and non-accreting pre-main sequence 
stars show any systematic difference in the abundance pattern. 
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and to investigate possible dependences of the coronal abun- 
dances on the activity level. 

The paper is organized as follows. Section |2] summarizes 
relevant information about the TMC and the XEST observa- 
tions. In Sect. |3]we describe the strategy employed to select the 
sample and we report the spectral fitting results, whose discus- 
sion is presented in Sect.|4] Finally, we draw our conclusions in 
Sect. |5] Appendix |A] details the simulations performed for this 
work. 

2. The TMC and the XEST observations 

The TMC is a cloud complex at 140 pc from the Sun, which 
contains about 400 known low-mass members; they are spread 
over a relatively large portion of the sky (~ 100 square degrees) 
and their ages cover a range between 10^ and 3 x 10^ years, 
with most stars being typically 1-10 million years old. The 
XEST survey consists of 19 recent observations (fields XEST- 
02 to XEST-20), obtained in two separate periods (August- 
September 2004 and February-March 2005) and with expo- 
sures ranging from 27 to 35 ks, to which 8 more previous fieldll] 
from the XMM-Newton archive were added, having longer ex- 
posure times in the range 43-131 ks. In total, the surveyed area 
covers only ~ 5 square degrees of the Taurus Molecular Cloud; 
however, since the observed fields were chosen so as to include 
the densest regions of the TMC with the largest concentrations 
of stars, the total number of detections associated to Taurus 
members in the EPIC images is 159. Including non-resolved 
companions that might be contributing to the X-ray emission, 
the XEST survey likely detected around 210 members, that is 
about a half of the entire known stellar population in Taurus. 
This sample contains mainly classical and weak-lined T Tauri 
stars, but also protostars, brown dwarfs and a few Herbig Ae/Be 
stars and other members with uncertain classification. 

The present study is based on X-ray spectra taken with the 
most sensiti ve instrument onbo ard XMM-Newton, the EPIC 
PN camera (IStruderetalJl200lh . This non-dispersive CCD de- 
tector has an effective area of ~ 1200cm^ at 1.5 keV and pro- 
vides moderate spectral resolution R = E/AE ~ 5 - 50 in the 
range 0. 1 - lOkeV. Bright sources in the field XEST-26 (around 
SU Aur), which lacks PN data, were excluded from the anal- 
ysis to ensure homogeneity of the results. For the XEST ob- 
servations, the PN was operated in full window mode and the 
medium filter was applied in all observations, with the excep- 
tion of the field XEST-27 (around ^ Per) where the thick filter 
was used. 



We refer to the paper bv Giidel et al.l ( 12007b for more details 
about all the observations of the XEST project, which were 
processed using SAS v6. 1, and the data reduction steps. 

3. Data analysis 

3.1. Preliminary simulations and choice of the sample 

A study of the chemical composition and of possible effects 
due to activity or accretion processes requires a large sam- 



' The field L1495 around V410 Tau was observed twice; the two 
observations are labelled XEST-23 and XEST-24. 



pie of stars and accurate and precise abundance measurements. 
Accurate estimates of abundances can be derived with the anal- 
ysis of high-resolution spectra, where individual emission lines 
of several chemical elements are well resolved. The XEST sur- 
vey has provided useful RGS spectra for nine TMC members, 
which have been analysed by iTelleschi et al.l (l2007h . 

In order to consider a larger sample of T Tauri stars, we 
based our study on the medium resolution spectra taken with 
the PN camera, which is ~ 10 times more sensitive than the 
RGS spectrometers and, unlike the latter, can also provide data 
for sources located off-axis in the field-of-view. The use of PN 
data required a preliminary feasibility study to understand the 
reliability of the results obtained from their analysis. In fact, the 
blending of the lines of many elements (neon, iron, nickel, oxy- 
gen) and the method used to analyze these spectra (i.e. global 
fitting with a plasma emission model with few isothermal com- 
ponents) might lead to systematic uncertainties in the abun- 
dances derived for some of the elements. In particular, we ad- 
dressed two issues: (i) Is the global fitting of PN spectra able 
to recover the true values of elemental abundances, consider- 
ing sources in a variety of physical conditions? (ii) What min- 
imum number of detected source counts are required to obtain 
measurements with relatively low abundance uncertainties, say 
within a factor of 2? 

These two issues are clearly not unrelated, since the ability 
of deriving an accurate and precise abundance value for a given 
element depends on a) whether the lines of that element have a 
good signal, and hence on the element abundance, on the ther- 
mal structure of the plasma and also on the amount of interstel- 
lar absorption, b) the statistics of the entire spectrum, which 
must allow a good determination of the source thermal prop- 
erties. We also point out that facing these issues requires the 
most accurate knowledge of the instrumental response function 
allowed by the instrument calibration. 

To answer the above questions, we performed a large num- 
ber of simulations exploring a variety of source conditions. At 
the same time, this procedure allowed us to put a rough thresh- 
old on the total number of source counts, below which the 
fitting parameters begin to be poorly constrained, and hence 
to define the sample of stars t o be studied. The ap proach is 
analogous to that employed by Magg io et al.l (120071) for their 
abundance study on a sample of stars in Orion, observed with 
Chandra/AClS. 

In brief, our simulations, whose details are reported in 
AppendixlAl showed that the abundances of most elements can 
be rather accurately derived from fitting of PN spectra and set 
an indicative threshold of ~ 5000 total counts in order to have 
useful PN data for abundance studies. Uncertainties at 90% 
level are usually better than a factor ~ 2.5 for spectra with 
~ 5000 counts, and better than a factor ~ 2 for spectra with 
~ 10000 counts. The iron abundance is the best constrained 
in most simulations, although some caveats must be noted in 
a few cases (see Appendix lAll. Instead, the abundances of Ar 
and Ca are generally very poorly constrained, as well as the O 
abundance in case of high absorption (A^e > 10^' cm"^). 

Therefore, we selected all the TMC members with more 
than 5000 counts in the PN spectra (0.3 - lOkeV), integrated 
over time intervals not affected by high background levels. We 
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Table 1. Sample of TMC members selected for the abundance 
study. "XEST id" is the identification label in the XEST cat- 
alog, the first two digits referring to the XEST field. For each 
source, both the total exposure time of the relevant observa- 
tion and the screened time (GTI) selected for this study are 
reported; "PN cts" are the source counts during the GTIs. "W" 
and "C" stand for WTTS and CTTS, respectively. TTS types 
and spectral types are taken from lGiidel et al.l (l2007l) and refer- 
ences therein, except where noted. 



Name 


XEST id 


Total 


GTI 


PN 


TTS 


Spec 






Exp. 




cts 


type 


type 






(ks) 


(ks) 








Anon 1 


20-005 


26.8 


26.8 


4520 


W 


MO 


CoKu Tau 3 


12-059 


26.4 


25.7 


10454 


W 


Ml 


DNTau 


12-040 


26.4 


25.7 


5169 


C 


MO 


IT Tau 


18-030 


26.8 


23.8 


8605 


C 


K2 


L1551-51 


22-089 


48.9 


45.4 


4506 


W 


K7 


V773 Tau 


20-042 


26.8 


26.8 


29971 


W 


K2 


V826 Tau 


22-100 


48.9 


45.4 


16950 


W 


K7 


V830 Tau " 


04-016 


27.6 


26.1 


8765 


w 


K7 


DH Tau " 


15-040 


27.4 


23.3 


14780 


c 


Ml 


HD 283572 


21-039 


41.1 


33.9 


79424 


w 


G5 


HP Tau/G2 


08-051 


33.7 


31.2 


15033 


w 


GO 


Hubble 4 


23-056 


62.6 


51.9 


16394 


w 


K7 


KPNO-Tau 15 


08-043 


33.7 


31.2 


6017 


w 


M2.75 


RYTau 


21-038 


41.1 


33.9 


5685 


c 


Kl 


V410 Tau 


24-028 


36.9 


27.4 


26740 


w 


K4 


V819 Tau 


23-074 


62.6 


51.9 


5524 


w 


K7 


V892 Tau 


23-047 


62.6 


51.9 


10886 


c 


K2 


BP Tau * 


28-100 


116.4 


88.0 


17770 


c 


K7 


TTau" 


01-045 


67.2 


37.5 


29960 


c 


KO 


RXJ0422. 1 + 1934 


01-054 


67.2 


37.5 


11296 


w 





" The light curve suggests that the emission is affected by the 

final stage of a flare decay. 
* A flare occurred during the observation is excluded from the 

GTIs. 

^ Secondary of the binary system. 



also added to the sample the two stars Anon 1 and LI 55 1-51 
with ~ 4500 counts. The sample defined in this way is listed 
in Table [1] and consists of 20 TMC members, 13 of which 
are weak-lined T Tauri stars, and the remaining 7 are accret- 
ing PMS stars (1 Herbig Be star and 6 CTTSs). Note that the 
Herbig Be star in this sample, V892 Tauri, has a close compan- 
ion (0.05" separat ion) with an estimated mass of 1 .5 - 2 Mq 



jSmith et alj 120051) . which therefore is likely to be the main 
source of X-ray emission, rather than the Herbig Be star itself. 
It is also likely that su ch a companion is an a ccreting star, based 



on previous studies ("Puch ene et all 1 19991: lK5nig et all 12001 



Itiartigan & Kenyon 2003.) showing that mixed binary systems, 
i.e. systems that include both accreting and non-accreting stel- 
lar compo nents, are v ery uncommon. Using the evolutionary 
models by Siess et al] ( 2000) and age equal to that of the pri- 
mary star (~ 3 Myr), the spectral type of this companion is K3 
for M = 1 .5 Mo, K2 for M = 1 .7 Mq and Kl for M = 2 Mq. In 
the following we will assume a spectral type of K2. 



3.2. Spectral fittings of TMC members 

The PN counts of the selected TMC members were extracted 
from circular regions around the sources, having typical radii 
of ~ 40" - 50". For each source, the background events were 
extracted from an annulus around the source; contamination 
fror n close X-ray sour ces within the annulus was also avoided 
(see' Giidel et al .1 12007 1 for more details about the calculation of 
the source and background extraction regions). The standard 
tasks of SAS were used to build the instrumental responses for 
all sources. In Table[T]we report both the total exposure times 
of the fields where the selected targets lie and the Good Time 
Intervals obtained by excluding from the present analysis those 
time intervals affected by high background emission (due to 
solar events), in order to increase the S /N ratios of the stud- 
ied spectra at high energies; in the case of BP Tau, a ~ 30 ks 
interval when the source exhibited a strong flare was also ex- 
cluded, so as to avoid possible effects of time variation of 
the abundances, as suggested by a number of previous studie s 
( Tsuboi et al J 1998HGudel etalJl999l: iFavata & Schmittll999l) . 

All PN spectra were fitted with the software XSPEC, lim- 
iting the fitting of each source to those energy bins where the 
background spectrum does not affect significantly the source 
emission (usually 0.3 - 7keV, in very few cases we excluded 
the energy bins at £ > 5 keV). We used a two-component ther- 
mal, optically thin plasma based on the APEC emissivity code, 
plus interstellar absorption, which is the same model employed 
for the simulations. The model has 13 free parameters: the hy- 
drogen column density A^h, the temperatures Ti and T2 and 
the emission measures EMi and EM2 of the two components, 
and the abundances of O, Ne, Fe, Mg, Si, S, Ca and Ar. For 
each source, the A^h value, the two temperatures and the ratio 
of the emission measures of the best-fitting model are reported 
in Table 121 together with the unabsorbed X-ray luminosity de- 
rived from the model, the average temperature and the surface 
X-ray flux derived from Lx and the stellar radius (the latter is 
taken from the Tables in lGiidel et alJl2007h . Average tempera- 
tures are weighted by the emission measure and are useful to 
compare the global thermal properties of the emitting plasma in 
different stars; they are indicative values and are reported with- 
out errors in the table. The surface fluxes are also reported with- 
out errors, owing to uncertainties in the adopted stellar radii. In 
Fig. [T]we plot the observed spectra with their relevant model 
spectra, for some example cases. 

Figure |2] shows the properties of the X-ray emitting plasma 
of the studied sample. The distributions of temperatures and 
EM ratios are very similar in the subsamples of WTTSs and 
CTTSs; the latter seems to have, on the whole, average temper- 
atures slightly greater with respect to the WTTS (the medians 
of the distributions of T^y are 13.5 MK and 19.5 MK, respec- 
tively), although a Kolmogorov-Smirnov test gives 7% proba- 
bility that the distributions for WTTSs and CTTSs are not dis- 
tinguishable. 

In Fig. |3]we plot the derived abundances (in solar units of 
Anders & Grevess3ll989 ) of the TMC members, where the el- 
ements are ordered by increasing FIP. The 20 values of abun- 
dances are rather widely distributed. The medians of their dis- 
tributions show a pattern with a minimum around Fe and Si, 
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Table 2. Parameters of the best-fit models for the 20 TMC members, with Icr errors. The unabsorbed X-ray luminosity Lx and 
the surface flux Fx are calculated in the 0.3 - lOkeV band. T.^^ is defined as T^^ = {EMiTi + EM2T2) / (EMi + EM2). 



Name 


Nii 


Ti 


T2 




Xv2 (d.o.f.) 






Fx 




(10-' cm--) 


(MK) 


(MK) 






(10^" erg s"') 


(MK) 


(lO'^ergs"' cm"-) 


Anon 1 


4.9 ±0.7 




29+f 


0.12 ±0.06 


0.96(141) 




7.4 


10.0 


CoKu Tau 3 


4.6 ± 0.4 


4.3 ± 6.7 


18.6 ± 1.2 


1.2 ±0.4 


0.88 (284) 


7.4+2i| 


12.2 


20.6 


DNTau 


0.44 ± 0.20 


8.7 ±0.8 


24 ±3 


1.6 ±0.8 


1.18 (142) 


1 10+"'^ 


18.4 


3.4 


IT Tau 


7 0+'- 

-0 6 


nrl 


35 ±5 


4.4 ±3 


0.92 (256) 


6.6!tr 


30.7 


22.8 


L1551-51 


0.78 ± 0.23 


7.4 ±0.7 




0.6 ±0.3 


1.04(128) 


1 75+''-25 

■ -0.09 


12.3 


14.0 


V773 Tau 


1.58 ±0.10 


4 8+°-5 


24.3 ± 1.3 


2.0 ±0.3 


1.07(434) 


9 0+' ^ 
'■"-0.6 


17.9 


40.2 


V826 Tau 


1.03 ±0.11 


4 7+0.8 
^- ' -0.4 


17.9 ± 1.6 


0.89 ± 0.22 


1.16(304) 


5.1 ±0.3 


11.0 


21.6 


V830 Tau 


0.79 ± 0.20 


6.8!^;^ 


23 ±6 


0.9 ± 0.4 


1.17(237) 


5 7+'° 
'-0.4 


14.5 


28.7 


DHTau 


1.82 ±0.20 


8.9 ±0.5 


25 ±3 


1.1 ±0.4 


1.04(319) 


7 7+0.9 
'■'-0.4 


17.3 


37.7 


HD 283572 


0.55 ± 0.05 


8.5 ±0.2 


24.5 ± 1.1 


1.96 ±0.21 


1.02(599) 


12.4 ±0.3 


19.1 


30.7 


HP Tau/G2 


3.6 ± 0.3 


8.7 ±0.5 


22.8 ± 1.4 


1.9 ±0.5 


0.94 (329) 




18.0 


24.1 


Hubble 4 


2.96 ± 0.21 


4 7+0.6 
^■-'-1.5 


19.2 ± 1.7 


1.01 ±0.21 


1.14(341) 


6 2+'^-^ 
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Figure 1. Examples of PN spectra (points with error bars) and best-fit model spectra (histogram) of selected TMC members in 
the present study. The bottom of each plot shows the residuals in units of x- 



a trend already obtained from the analysis of high-resolution 
spectra for several main sequence, pre-main s equence and gi- 
ant s tars with h igh X -ray luminosi ty CSanz -Forcada et al. I2OO3I; 
Argir offi et all 120041: ISc elsi et all 120041: lArgirofli et al. 2005; 
Scelsi et alj uOOSt iTelleschi et alj l2005h . However, as men- 
tioned above, the Ca abundance was generally ill-constrained 
in the simulations, which showed that nearly solar values can 
be frequently obtained even when the input Ca abundance is 
significantly subsolar; hence, if we do not consider calcium, 
the pattern of the medians resembles an inverse FIP effect. 



It is important to note, however, that the medians of the 
abundance distributions express average properties of the stud- 
ied sample, while the abundances of individual stars might be 
correlated with other stellar parameters. Yet, comparison of 
Fig.[3]with the results of the simulations shows that the spread 
of the derived abundances for a given element is roughly of 
the same magnitude as the one found from the simulations (in 
comparing Fig|3]and Fig. |AT| consider that half of our sample 
stars have less than ~ 10000 counts). This implies that with 
the available data we can hardly distinguish between the abun- 
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T, (MK) T2 (MK) EM2/EM, T„, (MK) Lx (10^" erg s"') 

Figure 2. Temperature, emission measure and X-ray luminos- 
ity values for the WTTSs and the CTTSs of the studied sample. 
The values on the y axis have different units, for each quantity, 
as reported in the x-axis labels. 
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Figure 3. Abundances of the 20 TMC members for each ele- 
ment, ordered for increasing FIR The size of the circles gives a 
rough indication of the accuracy of the measurements: large 
circles for abundances derived from spectra with less than 
10000 counts, medium circles for 10000 < counts < 20000, 
and small circles for counts > 20 000. At the bottom of the plot, 
close to the label of each element, the number of stars with 
abundances lying outside the plot (i.e. lower than 0.01 solar) 
are indicated. The short thick segments mark the medians of 
the abundance distribut ions, while th e long thin segments mark 
the values assumed by IGiidel et al.l (12007.) to fit EPIC spectra 
of Taurus members (see Sect. 14.1b . 



dances of any given element in different stars of the present 
sample, and hence it is not meaningful to search for possible 
correlations with other parameters. This result also suggests 
that the derived sets of abundances of these TMC members iMa ggio et al. 2007 
could be drawn from a single set of abundances, i.e. the spread .2001; Audard et al.l 2003h . 



of the values for any given element is not due to diff'erences 
among the sources, but rather to statistical fluctuations. To 
check this idea, we proceeded as follows: we simulated 10000 
PN spectra starting from a two-temperature model which re- 
flects the average observed properties of the studied sample, 
that is A^H, Ti, T2, EMojEMi and elemental abundances equal 
to the medians of the relevant observed distributions (TableO. 
The total number of counts in the 10000 simulated spectra re- 
flects the distribution of counts observed in the real spectra of 
the 20 TMC members. Afterwards, we fitted these spectra with 
a 2-T model with the same (13) free parameters as before. 

The result of this test is reported in Fig. H] which shows 
the distributions of the abundances (in the upper panel) and the 
abundances normalized to iron (in the lower panel) obtained 
from the real data, compared with the central 90% ranges of 
values obtained from the simulations, for each element. In both 
panels, the spread of the values derived from real data is com- 
parable to that indicated by the simulations, for all elements. 
Exception is the Ne/Fe abundance ratio, whose distribution is 
significantly larger in the case of the real data. This is better 
quantified by a K-S test which shows that for the case of the 
Ne/Fe ratio the two distributions are diff'erent with a signifi- 
cance of more than 3cr. The range of O/Fe ratios also appears 
to exceed the 90% bounds of the range of the simulations, but 
the difference between the two distributions is not statistically 
significant (K-S probability of undistinguishable distributions 
equal to 7%). The result about the Ne/Fe ratio is intriguing 
and deserves to be further investigated, because neon and iron 
are chemical elements having different properties related to 
their tendency of forming dust grains, with the volatile noble 
gas neon remaining in the gaseous phase while iron more eas- 
ily condensing in grains. Moreover, Ne is a high-FIP element, 
while Fe is a low-FIP element. In the next Sect. I4.2l we discuss 
the Ne/Fe ratio for the stars of the present sample together with 
some possibilities that could explain the above result. 

In Fig. |5] we separate the abundance distributions of the 
classical and the weak-lined T Tauri stars: the plot shows that 
there is no significant difference between the distributions for 
the two classes of sources, indicating that both the thermal 
properties (Fig. |2]l and chemical composition of the X-ray- 
emitting plasma are similar 



4. Discussion 

4.1. Medians of the abundance distributions 

A robust result of this analysis is the low iron abundances 
of the X-ray bright Taurus members with respect to the so- 
lar photospheric value (Fe ^ 0.2 solar), since the abundance 
of this element is generally well-determined and with rela- 
tively small uncertainties (as indicated by the simulations). 
This analysis thus confirms the results of other previous stud- 
ies on high-luminosity stars (pre-main sequence stars and 
RS CVn binary systems) which indicate a low coronal metal- 
licity, with respect t o the solar case, in s uch active stars (e.g. 



Drake et al. 2001 ; Huenemoerder et al 
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Table 3. Medians of the distributions of absorption, temperatures, emission measures and abundances derived from the PN 
spectra of the studied sample. Abundances are in solar units, A^h in 10^' cm Ti and T2 in MK. 
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Figure 4. Distributions of abundances {upper panel) and abun- 
dances relative to iron (lower panel) found for the 20 TMC 
members, compared with the results of simulations (see text). 
The boxes contain the central 90% of the abundance values ob- 
tained from the simulations, the median being indicated by the 
thick black segment within each box. The thin, longer segments 
mark the medians of the abundance distributions derived from 
real data; they are practically coincident with the median val- 
ues of the simulations. Numbers at the bottom of the plot are as 
in Fig. [3] 



In Fig.[3]we also r eport the fixed values of abundances used 
bv lGiidel et alj (120071) for fitting the EPIC spectra of more than 
a hundred Taurus members, with the aim of deriving absorp- 
tion column densities, characteristic plasma temperatures and 
X-ray luminosities. The values adopted by these authors (see 
also Appendix lAb were chosen as representative of the coronal 
abundances in active stars. Our detailed analysis confirms the 
validity of this assumption, because the medians of the abun- 
dance distributions for all elements are very similar to the pre- 
viously adopted values, with the possible exception of calcium. 

In Sect. 13.21 we showed that the medians of the abun- 
dance distributions follow a pattern vs. FIP similar to that 
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Figure 5. Scatter plot of abundances obtained for the classical 
(open circles) and weak-lined (filled circles) T Tauri stars in the 
studied sample. Numbers at the bottom of the plot are as in Fig. 
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Scelsi et al.l l2005t iTelleschi et al. 2005); given the large uncer- 
tainties of the Ca abundances, our analysis might also suggest a 
systematic pattern of increasing abundances for increasing FIP. 
However, such a pattern is obtained by adopting solar photo- 
spheric abundances as reference, while any FlP-related effect 
should be studied by comparison with the stellar photospheric 
abundances; unfortunately, they are not available for any of the 
stars in our sample. Measuring abundances in stellar photo- 
spheres is a difficult task, due to the limited quality of stel- 
lar optical spectra, to uncertainties in stellar parameters (grav- 
ity, effective temperature, macro-turbulence) and atmospheric 
models (LTE approximations), to line blending (especially for 
cool M-type stars) and rotational broadening effects which is 
typically the case of active young stars; furthermore, the ef- 
fects of possible stellar spots and of the chromospheric activity 
on abundance determinations are not well understood, and bi- 
narity requires the non-trivial procedure of untangling compos- 
ite spectra. An assumption usually adopted is that stars in the 
solar neighbourhood have photospheric abundances similar to 
the solar ones, because they belong to the same galactic pop- 
ulation; moreover, we expect that stars younger than the Sun 
have on average higher metallicities because they have formed 
from material enriched in high-Z elements by the remnants of 
an older stellar population. However, extensive surveys show 
that there is a large scatter in metallicity among stars of similar 
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age dNordstrom et al.ll2004 . ISanz-Forcada et all ( |2004 stud- 
ied the coronal abundances of four stars with accurately known 
photospheric abundances and showed indeed that the pattern 
of coronal-to-photospheric abundance ratios vs. FIP can be 
different if stellar o r solar photospheric abundances are used. 
The recent study by 'Maggio et alJ ( 2007 ) on the X-ray bright- 



est stars in Orion also indicates that a clear trend vs. FIP is 
no longer observed when the photospheric abundances of the 
Orion stars or the composition of the gaseous phase in the neb- 
ula are considered instead of the solar ones. On the other hand, 
ISanz-Forcada et al.l (l2003h derived coronal abundances that are 
not photospheric for AB Dor, which h as well measured photo- 
spheric abundances, and iTelleschi et al. (2005) found the FIP 
and inverse FIP patterns in six solar analogs for which sev- 
eral photospheric abundances are rather clearly known. It is 
important to dedicate more observing time to collecting optical 
spectra in order to measure stellar photosperic abundances and 
hence to address in a proper manner the issue of any FlP-related 
difference between the photospheric and coronal composition. 

4.2. Ne/Fe abundance ratio 

We now return to the spread of abundances found for the sam- 
ple stars, for each element. The only quantity showing a spread 
of the values larger than expected from the simulations is the 
Ne/Fe ratio (Fig.|4|i, which varies between ~ 2 and ~ 9, with a 
median around Ne/Fe~ 5. The ratio between the abundance of 
a given element and that of iron is generally better constrained 
than absolute abundance values, because systematic errors in 
the elemental and iron abundances, due e.g. to a wrong esti- 
mate of the continum emission, usually go in the same direction 
and can compensate when calculating their ratio. The result we 
obtained deserves to be investigated for several reasons. 

A plausible assumption is that the photospheric abundances 
of iron and neon are about the same in all the stars studied here, 
because they belong to the same region and hence they were 
born from the same cloud material; then, the Ne/Fe ratio in 
their coronal plasma could vary depending on the stellar activ- 
ity level if FlP-related effects are really present in stellar coro 



nae, as s uggested by previous studies (IBrinkman et al.l 12001 



Audard et al . 2003). The X-ray luminosities of the sample stars 
are in the range 10"*° - 10^' erg s theref ore the study is biased 
towards active TMC members. However. I Audard et al.l ( 2003 ) 
found that a transition from an inverse FIP effect to a flat pat- 
tern with decreasing coronal activity is visible in a sample of 
RS CVn binary systems spanning about the same X-ray lumi- 
nosity range as that of the present study. In particular, they re- 
port an anticorrelation between Fe/O vs. the average coronal 
temperature T^y of five active RS CVn-like systems as well as 
a constant value of the Ne/O ratio, which translates into a cor- 
relation between N e/Fe and ray . When adding a sample of solar 
analogs studied bv lGiidel et al. I (l2002l) . reaching lower activity 
levels, such a correlation becomes much more evident and a 
trend from inverse to normal FIP effect is seen. 

We plot in Fig. |6]the Ne/Fe ratio for our sample of TMC 
stars as a function of various activity indicators (X-ray lumi- 
nosity, Lx, average temperature, T.^y, and surface X-ray flux. 



Fx), in order to search for possible transitions from a FIP ef- 
fect to an inverse FIP effect with increasing activity. However, 
no significant correlation of Ne/Fe with any of the considered 
parameters is found from the available data, which therefore do 
not support the existence of different FlP-related effects in the 
coronae of the studied stars. 

However, we also note an apparent drop of the Ne/Fe ra- 
tios for stars with average temperature T^y > 15 MK. Based 
on the analysis of RG S spectra of a few stars in the XEST 
project, ITelleschi et al. 1 12007) also found a bimodality in the 
Fe abundance and Ne/Fe ratio for their sample stars, although 
this behaviour was observed to be linked to the spectral type. 
In particular, the above authors found larger Fe abundances and 
lower Ne/Fe ratios in three G-type stars with respect to six K- 
and M-type stars. To check their result, we divided our sample 
into two groups, i.e. stars with spectral type earlier and later 
than K2 (7 sources belong to the first group, while 12 stars 
have spectral type later than K2, the spectral type of the star 
RXJ 0422. 1-1- 1934 being not available), and we observe in Fig. 
|6]that most of the stars cooler than 15 MK have indeed spec- 
tral type later than K2. The K-S test gives 0.15% probabil- 
ity that the distributions of T-^y for the two groups are drawn 
from the same parent distribution. We noted that this difference 
derives from earlier-type stars having typically higher T2 and 
EM2IEMX values. 

Figure|7]shows the cumulative distributions of the iron and 
neon abundances and of the Ne/Fe ratio for the stars of the 
two groups: there is a systematic difference in the neon abun- 
dance against the spectral type, and an even more clear evi- 
dence for lower Fe abundances in K- and M-type stars which 
implies larger Ne/Fe ratios than in earlier-type stars. The prob- 
ability that the Fe distributions for the two groups are not 
distinguishable is 0.6%; for the cases of the Ne/Fe and Ne 
distributions, the probab ilities are 8% and 17%, respectively. 
Franciosini et al. ( 200d) studied the G-type Taurus member 
SU Aur (r^v ~ 20 MK) and found Fe ~ 0.6 times the solar 
value and Ne abundance about equal to the solar one, with in- 
dependent analyses of MOS and RGS data. If we consider also 
the abundances obtained for SU Aurigae in this picture, the re- 
sult for Fe and for Ne/Fe is statistically more significant; the 
three K-S test probabilities mentioned above are 0.4% (Fe), 
4.5% (Ne/Fe) and 29% (Ne). In this calculation we have as- 
sumed that the spectral type of the companion of V928 Tau 
(the probable source of X-rays) is K2; we have verified that 
the above conclusions do not change assuming for this star a 
spectral type later than K2. 

A possible dependence of Fe and Ne/Fe on the spectral 
type has not been cl aimed in the past, before the report by 



Telleschi et alJ (l2007h . and high-quality spectra of a larger sam- 



ple of stars are needed to confirm it. Since the stars in our 
sample with spectral type later than K2 have, on average, 
cooler coronae, this result may also reflect a dependence on 
Tav (although this is not evident in the Telleschi et al. sam- 
ple). However, for stars with the highest T^^ (in the range 
10 - 30 MK), such a dependence wou ld be in the opposit e 
direction with respect to that found by lAudard et al.l (l2003h . 
who report increasing Ne/Fe with increasing T^y in the range 
~4 - 15MK. 
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Figure 7. Ne and Fe abundance distributions and Ne/Fe ratios for stars with spectral type earlier and later than K2 (three upper 
panels); the circle size decreases with increasing number of total counts in the PN spectrum. In the three lower panels the relevant 
cumulative distributions are shown. 



In the light of the above result, we run simulations anal- 
ogous to those described in Sect. 13.21 but taking into account 
a possible bimodality of the iron abundance with the spectral 
type, i.e. we assumed Fe= 0.24 and Fe= 0.15 for stars with 
spectral type earlier and later than K2, respectively. Figure [8] 
analogous to the lower panel of Fig. ID shows that the spread 
of Ne/Fe in real data can be explained in this way (K-S test 
probabiUty = 32%). 

Another issue of potential importance is that Fe forms 
grains much more readily than Ne: as explained in Sect. [T] if 
iron settles in dust grains in the disk of CTTSs and the dust and 
gas components are separated for some reason. X-ray emis- 
sion from accreting stars may reveal a metal-depleted plasma. 
We note that such a result may be expected if accretion-driven 
emission dominates over a coronal component, since in this 
case the abundances we are measuring are essentially those of 
the accreted plasma. 

The X-ray emitting plasmas of the stars in our sample have 
rather high average temperature (> 10 MK), and similar prop- 
erties in both the CTTSs and the WTTSs, indicating that we are 
probably seeing a dominant contribution from a coronal emis- 
sion. Figure |9] shows the Ne/Fe ratio of our sample stars as 
a function of two stellar parameters linked to the presence of 
circumstellar material, i.e. A^h and the accretion rate. For com- 
parison, we report also the ratios for the two classical T Tauri 
stars TW Hvdrae dRobrade & Schmittll2006t) and MP Muscae 
dArgiroffi et al.ll2007h : for the CTTS BP Tauri, previously stud- 
ied by other authors, we note tha t our Ne/Fe ratio is in ver y 
good agreement with that found bv lRobrade & SchmittI (120061) . 



As expected, we do not find evidence of any relation between 
Ne/Fe and these two parameters; moreover, as shown in Fig.|9] 
high Ne/Fe ratios may be present both in CTTSs and in WTTSs 
and are therefore not a peculiar characte ristic of accreting stars . 
This is consistent with the findings of ITelleschi et al.l (l2007h 
that high Ne/Fe ratios are not present in all accreting stars. 

4.3. Ne/0 ratio 

The Ne/O ratio shares some similarities with the Ne/Fe 
one, in that oxygen can easily combine with dust grains; 



Drake & Testa ( 20051) analyzed high-resolution Chandra spec- 
tra of a variety of stars (dwarfs, giants, multiple systems, 
subgiants) and found that the Ne/O ratio is remarkably con- 
stant in stellar coronae (Ne/O ~ 0.4, the solar value by 
Anders & Grevessdll989l being ~ 0.14). This is true also for 
the classical T Tauri stars BP Tauri and TWA 5, but not for 
the CTTS T W Hydrae, whose h igher Ne/O value (~ 1) was in- 
terpreted bv lDrake et alj (l2005h in terms of depletion of grain- 
forming elements in the X-ray emitting gas which is accret- 
ing onto the star: such a depletion is supposed to be more pro- 
nounced in the case of the more evolved and older circumstel- 
lar disk of TW Hydrae, where coagulation of grains into larger 
bodies orbiting the star could have happened; in the younger 
CTTS BP Tauri, where part of the emission is accretion-driven, 
the dust particles could be still small and can release, after 
heating, the elements l ocked into the grains. However, very 
recently lArgiroffi et"ai] ( [2007.) obtained Ne/O ~ 0.5 for the 
CTTS MP Muscae, which is even older than TW Hydrae, but 



10 



L. Scelsi et al.: Coronal abundances in Taurus-Auriga X-ray briglit stars 




" 10.0 



Fx (10^ erg s"' cm"^) 



Figure 6. Ne/Fe ratios of all sample stars plotted as a function 
of three activity indicators: X-ray luminosity (upper panel), av- 
erage temperature [middle panel) and surface X-ray flux [lower 
panel). In the plot of Ne/Fe vs. Tav, different spectral types are 
shown (see text): stars with spectral type later than K2 are plot- 
ted as filled circles, while open circles refer to spectral type ear- 
lier than K2 (an open square superimposed to the circle marks 
the Herbig Be star V892 Tau); the cross represents the WTTS 
RXJ0422.1-H1934, whose spectral type is not available. 



still shows evidence for a substantial contribution to the X- 
ray emission from accretion shocks in the high plasma den- 
sity derived from the analysis of the O vii triplet; the average 
plasma temperature of this star is relatively high and hence the 
derived abundances are likely afifected by the coronal emis- 
sion. What is also emerging from these abundance studies is 
the need for a criterion to understand which abundances (coro- 
nal or accretion-related) we are measuring in these stars, per- 
haps also with the help of modeling. For the stars we studied 
here, the Ne/O ratio is affected by large uncertainties in those 
sources with high absorption, because of the difficulty of con- 
straining the O abundance in this case. If we consider the 12 
stars with A^h < 3 x 10^' cm"^, we find that the Ne/O ratio is 
in the range ~ 0.2 - 0.6, with the values equally spread around 
~ 0.4 (Fig. [Tol l. In particular, this is true also for the classi- 
cal T Tauri stars DN Tauri, DH Tauri, T Tauri and BP TaurjJ 
among the 12 we are considering. They are rather young, with 
ages in the range ~ 1 - 2.7 Myr, and their high average plasma 
temperatures imply a large coronal contribution, as already dis- 
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^ The value Ne/O found in this study for BP Tauri i s in good agree- 
ment with that obtained bv lRobrade & SchmittI JioO^ . 



Figure 8. Distributions of abundances relative to iron for the 20 
TMC members compared with the results of simulations taking 
into account a possible bimodality of the iron abundance (see 
text). Symbols and lines are as in Fig.|4] 



cussed. Therefore, our result does not go against the scenario 
mentioned above. 



5. Summary and conclusions 

The thermal properties and chemical composition of the X-ray 
emitting plasma in a sample of bright Taurus members have 
been studied in this work with the aim of investigating pos- 
sible systematic differences between classical and weak-lined 
T Tauri stars, and of understanding whether the pattern of abun- 
dances in such stars is related to some stellar or environmen- 
tal characteristics. We used data from the XEST survey; since 
there are very few useful RGS spectra of Taurus members, 
we used the PN data in order to enlarge as much as possible 
the sample to be analyzed. Preliminary extensive simulations 
showed that PN spectra with ~ 5000 counts or more in the 
0.3 - lOkeV band can provide reliable estimates of elemental 
abundances and 68% statistical errors usually within a factor 
~ 2 for most of the elements; exceptions may be O, Ca and 
Ar in some cases, depending on the statistics and source ab- 
sorption. The sample of stars considered for this work there- 
fore consists of 20 TMC members: 13 weak-lined T Tauri stars 
and 7 accreting stars, having at least 4500 total counts in the 
spectrum. Their medium-resolution spectra were globally fitted 
with an absorbed two-temperature thermal model with variable 
abundances of O, Ne, Mg, Si, S, Ar, Ca and Fe. 

The temperatures of the two components were found in the 
ranges ~ 5 - lOMK and ~ 20 - 40 MK, and the hot- to-cool 
emission measure ratios range between ~ 0.8 and 4. Average 
temperatures are in the range ~ 8 - 30 MK, with the classical 
T Tauri stars being slightly hotter than the WTTSs. 

The spread of the 20 abundance values for each element 
(Fig. O is more or less compatible with the errors obtained 
from the simulations; a posteriori tests, also based on simu- 
lations, have shown that our results are consistent with these 
TMC stars sharing a unique set of abundances, with the possi- 
ble exception of iron. 
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Figure 10. Ne/O ratios as a function of the photoelectric ab- 
sorption, for stars with A^h < 3 x 10^' cm"^. Errors are Icr. 
Filled and open circles mark WTTSs and CTTSs, respectively. 
For comparison, the positions of the two CTTSs TW Hydrae 
and MP Muscae are also shown (open circles without error 
ba rs). The dotted hne m arks the average value of 0.41 found 
bv Drake & Testal (120051) for stars of various activity levels. 



Figure 9. Ne/Fe ratios as a function of the photoelectric absorp- 
tion (upper panel) and accretion rate (lower panel). Errors are 
Icr. Filled and open circles mark WTTSs and CTTSs respec- 
tively. In both plots, the positions of the two CTTS TW Hydrae 
and MP Muscae are shown as open circles without error bars. 
In the lower panel, only stars with measured accretion rates are 
shown; all WTTSs have upper limits indicated by left arrows. 
Accretion rates for the TMC members, for TW Hydrae and fo r 
MP Muscae are ta ken from the tables in iGiidel et al. ( 2007 ). 
from lAlencar & B atalha (2002,) and from Sacco G. G. (private 
comunication), respectively. 



The main conclusions from the present work are summa- 
rized as follows. 

Hot coronal plasma is responsible for most of the X-ray 
emission from both CTTSs and WTTSs studied here. The 
coronal plasma has similar properties in the two classes of 
T Tauri stars, in particular similar characteristic temperatures 
and chemical composition. Late-type stars therefore develop 
hot coronae early in their life (within ~ 1 Myr), whose proper- 
ties change little during the following pre-main sequence evo- 
lution, at least up to an age of ~ 10 Myr The reason why the 
coronal emission appears suppressed in some CTTSs, such as 
TW Hydrae, and not in most of the others is to be clarified with 
future work. 

The available data do not allow us to evaluate quantitatively 
the presence, in the CTTSs considered here, of a cool X-ray 
emission excess that might be present due to shock-heating. 

The iron abundance is significantly subsolar (~ 0.1 - 0.5 
times the solar photospheric value) in all the sample stars. Since 
the abundance of this element is the best-constrained one, this 



result is rather robust and confirms previous finding for PMS 
and active stars. 

We found indication of coronae in G-type and early K-type 
stars hotter than in stars with later spectral type. We also find 
strong support for a bi-modality of the Fe abundance (and the 
Ne/Fe ratio), depending on the spectral type (or perhaps on the 
average temperature), with higher Fe values and lower Ne/Fe 
ratios for G-type and early K-type stars compared to stars 
of later spectral type. The latter res ult supports the findi ngs 
from high-resolution spectroscopy bv lTelleschi et al.l (l2007h . A 
larger sample is needed, however, to find out if this trend is 
merely the result of the limited statistics presently available. 

The analysis of the abundance distributions of this sam- 
ple of bright TMC members did not reveal any significant dif- 
ference in the abundances related to the presence of accre- 
tion/circumstellar material. Ne/Fe abundance ratios between 
2 and 9 times the solar value are found both for the classi- 
cal andJhe_weak-linedT Tauri stars, confirming the findings 
by iTelleschi et al.l (l2007b based on a smaller sample of high- 
resolution X-ray spectra of T Tauri stars. Regarding the Ne/O 
ratio, we do not observe in the classical T Tauri stars a be- 
haviour analogous to that of TW Hydrae, all values being be- 
low ~ 0.6. 
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Appendix A: Simulations 

In this section we describe the simulations we run before fit- 
ting the real data in order to assess the performances of the PN 
regarding abundance determinations. This part of the investiga- 
tion can be summarized as follows: 



1. 



3. 



Fourteen input coronal models were chosen, each consist- 
ing of an EM distribution and a set of elemental abun- 
dances, plus interstellar absorption. 

Three different numbers of source counts were considered: 
-5 000, 10000 and 20000 counts in the 0.3 - lOkeV 
band of the PN spectrum. This choice is motivated by the 
fact that most of the brightest TMC members have total 
PN counts within the above range; moreover, previous fit- 
tings of Taurus members with less than 5000 counts in 
the PN, using 2-T models with a fixed set of abundances 
dGiidel et al.ll2007b already provided good descriptions of 
the PN spectra, a result which implies no need to treat the 
abundances as free parameters. 

For each simulated observation, i.e. for each coronal model 
at each number of total counts, 1000 PN spectra were simu- 
lated (using the standard PN response function and assum- 
ing no background contamination) which have a Poissonian 
distribution of the total number of counts. 
For each simulated observation, the 1000 PN spectra were 
fitted with an automatic routine, using a 2-T coronal model 
with the following starting values for the fitting parameters: 
abundances and A^h equal to the values of the input coronal 
model, Ti and T2 given by the results of interactive fittings 
of a few simulated spectra for each of the fourteen input 
models. 



Moreover, for the input models 1 and 2 (see Table lA.ll 
below) the fitting of simulated spectra was repeated for each 
of the three values of the total counts, starting from different 
initial values of A^h and temperatures, and initial abundances 
equal to the solar values: for each of these 6 simulations, 300 
spectra were automatically fitted by performing an extensive 
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exploration of the Nu, Ti and T2 space so as to simulate the ap- 
proach to real data. The results obtained in this way are practi- 
cally undistinguishable from those derived with the procedure 
described in item #4 above. For the input models 9 and 10, non- 
solar abundances as starting values for the fitting of simulated 
spectra were also used. Again, the results in these cases are 
very similar to those obtained with the procedure in item #4. 

Note that using the same instrumental response both to sim- 
ulate and to fit the spectra implies that this investigation leads 
to estimate the errors in the determination of the plasma pa- 
rameters due to the photon counting statistics and to the instru- 
ment spectral resolution. The method cannot provide estimates 
of possible systematic errors due to uncertainties in the instru- 
mental calibration and in the atomic parameters used to calcu- 
late theoretical plasma emissivities. 

The input coronal models were chosen so as to explore 
characteristic absorption values and physical conditions of the 
coronal plasma (in terms of HMD and abundances), among 
those commonly found in the literature for active T Tauri stars. 
For the abundances, some extreme cases were also explored 
(see below). The emission measure distributions we have con- 
sidered for these simulations are defined in the temperature 
range 10''K< T < lO^Kas: 



EMD oc 



j-p 



io''K<r<rpeak 



(A.l) 



Based on EMD reconstructions with hig h-resolution spec- 
tra of i ntermediate- and high-activity st ars by Laming & Drake 



Argiroffi et all (l2005l) : iTelleschi et al 



(I1999I): lArg iroffi et al. ( 2 061: IScelsi et all j2005h: 



(120051); iGiidel et al 



(l2007h . we chose the values of a and p in the range 2-5 and 



1-2, respectively, while two values of Tpeak were explored 
(10^^ K and 10^ K). Four sets of abundances were adopted: (i) 
the abundance set used by lGiidel et al.l (l2007h to fit the spectra 
of known TMC members (hereaft er "TMC set"): with respec t 
to the solar photospheric values by I Anders & Grevesse (1 19891) . 
the abundances of this set are: C=0.45, N=0.788, 0=0.426, 
Ne=0.832, Mg=0.263, Al=0.5, Si=0.309, S=0.417, Ar=0.55, 
Ca=0.195, Fe=0.195, Ni=0.195; (ii) all a bundances equal to 
0.2 tim es the solar photospheric values of I Anders & Grevesse 
1989h : (iii) solar photospheric values for all abundances. 



except for iron {- 0.2 solar); (iv) the "TMC set", except for 
iron (= 0.5 solar). We finally considered two absorption levels: 
"low" absorption with A^h = 10^' cm"^, and "high" absorption 
with A^H - 10^^ cm"^. The fourteen input models are listed in 
Table O 

The two-temperature thermal models used to fit the simu- 
lated spectra are based on the APEC emissivity code and have 
the following 13 free parameters: the interstellar hydrogen col- 
umn density (A^h), the temperature and emission measure of 
each plasma component, and the abundances of O, Ne, Fe, Mg, 
Si, S, Ar, Ca. We fixed the abundances of C, N and Ni to the 
respective values in the TMC set, because they are generally 
not well-constrained from PN fitting; in fact, the C and N lines 
fall at low energies where the PN resolution is quite poor (and 
moreover the PN calibration is probably uncertain), while the 



Table A.l. Input mode ls for simulatio ns. "TMC" is the abun- 
dance set used by iGiidel et al.l (2007) (see text). "Low" and 



"high" absorption mean A^h 
respectively. 



lO^'cm"^ andA^H 



10^2 cm^l 



Model 


log Tpeak 


a 


P 


Abundances 


Absorption 


1 


7 


3 


2 


TMC 


low 


2 


7 


3 


2 


TMC 


high 


3 


7 


5 


2 


TMC 


low 


4 


7 


5 


2 


TMC 


high 


5 


7 


2 


1 


TMC 


low 


6 


7 


2 


1 


TMC 


high 


7 


7 


3 


2 


0.2 solar 


low 


8 


7 


3 


2 


0.2 solar 


high 


9 


7 


3 


2 


solar (Fe=0.2 solar) 


low 


10 


7 


3 


2 


solar (Fe=0.2 solar) 


high 


11 


7 


3 


2 


TMC (Fe=0.5 solar) 


low 


12 


7 


3 


2 


TMC (Fe=0.5 solar) 


high 


13 


6.7 


2 


1 


TMC 


low 


14 


6.7 


2 


1 


TMC 


high 



Ni lines are blended with the (usually) stronger iron and neon 
Unes. 

As representative examples of the results of all simulations, 
we report in Fig. lA.ll those obtained for the abundances in 
terms of medians of the distributions of the 1000 best-fit values, 
and error bars corresponding to the central 90% of the distribu- 
tions, for ten different cases illustrating the effects of differ- 
ent input models, statistics, absorptions and choice of starting 
abundance values for the fittings. 

The results indicated that, for all the cases analyzed here, 
the derived abundances are generally in good agreement with 
the values of the input models. The medians of the abundance 
distributions obtained from fitting simulated spectra are, for aU 
elements, compatible with the input value within 90% statisti- 
cal errors, and within the 68% errors in the great majority of the 
cases. Some systematic effects must be noted in a few cases: Ne 
results overestimated by a factor ~ 2 in models 7 and 11, while 
Fe and O are underestimated by about the same factor in model 
13. 

As to the magnitude of the statistical uncertainties, errors 
at 90% confidence level within a factor ~ 2.5 are obtained for 
most of the elements using PN spectra with 5000 counts, and 
within a factor ~ 2 with 10000 counts. Exceptions are oxygen, 
argon and calcium in some cases (see below). Iron is the best- 
determined element, constrained within a factor ~ 1 .6 when ab- 
sorption is relatively low, already with spectra containing 5000 
counts. For the other elements, 90% errors on the medians can 
be a factor ~ 1.4 up to 3 (spectra with 10000 counts), depend- 
ing on the model. 

The abundance of oxygen is very poorly constrained (even 
for the highest statistics) when the A^h is of the order of 
10^^ cm"^, as one could expect because of the strong attenu- 
ation of the spectrum at low energies (£ < 1 keV; recall that 
the strongest lines of the He-like and H-like ions of oxy- 
gen lie at £■ ~ 0.6 keV). Regarding argon, its abundance can 
be better constrained in the case of high absorption; this can 
be easily understood by noting that, for fixed number of total 
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Figure A.l. Medians (diamonds) and central 90% of the distributions of the abundances found from eight sample simulations, 
for different models and/or statistics. The input values of abundances are shown as squares. The input HMDs are equal in models 
1, 2, 7 and 9, with log Tpeak - 7. The EMD of model 3 is steeper, while the EMD of model 5 is flatter. The EMD of model 13 
peaks at log Tpeak - 6.7. The input models have low absorption, with the exception of model 2. The two bottom panels show the 
results obtained when the starting abundances for the fitting are different from the "input value" of the relevant model. 



counts, a larger fraction of all the spectral counts goes in the 
Ar lines (~ 3 keV) when the low-energy region of the spectrum 
is strongly suppressed by the interstellar absorption. Instead, 
although the strongest Ca lines fall at ~ 4.1keV, i.e. close to 



the Ar lines in the high-energy tail of the spectrum, the abun- 
dance of this element is generally ill-constrained, unless the 
abundance is of the order of the solar value, or higher, and the 
statistics are high. Since the PN effective area is about the same 
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at ~ 3 keV and ~ 4.1 keV, the difference with the case of Ar is 
likely due to the lower emissivity of the calcium Unes, com- 
pared to those of the Ar Unes. 

It is important to note also that the results of the simula- 
tions are practically the same, in terms both of magnitude of 
the errors and agreement between medians of the abundance 
distributions and input values, regardless of the starting abun- 
dance values used to fit the spectra. 

Finally, note that the simulations described in this section 
assumed no background contamination, which may not be neg- 
ligible at the highest energies. For the bright stars discussed in 
this paper, we verified that only in a couple of cases the back- 
ground spectrum dominates over the source spectrum in the en- 
ergy range £ > 5 keV, which was therefore excluded from the 
fitting. To investigate the effect of the exclusion of the Fe xxv 
6.7 keV complex on the ability to constrain the iron abundance 
and the other parameters of the model, we repeated the simula- 
tions for the models 1 and 2 by fitting the simulated spectra in 
the energy range 0.3 -5 keV. The results obtained in these cases 
are very similar to those obtained employing a larger spectral 
range. 



